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ABSTRACT 

We present a calculation of the atomic and low-ionisation emission line spectra of pho- 
toevaporating protoplanetary discs. Line luminosities and profiles are obtained from detailed 
photoionisation calculations of the disc and wind structures surrounding young active solar- 
type stars. The disc and wind density and velocity fields were obtained from the recently 
developed radiation-hydrodynamic models of Owen et al., that include stellar X-ray and EUV 
irradiation of protoplanetary discs at various stages of clearing, from primordial sources to 
inner hole sources of various hole sizes. 

Our models compare favourably with currently available observations, lending support 
to an X-ray driven photoevaporation model for disc dispersal. In particular, we find that X- 
rays drive a warm, predominantly neutral flow where the OI 6300A line can be produced by 
neutral hydrogen collisional excitation. Our models can, for the first time, provide a very good 
match to both luminosities and profiles of the low-velocity component of the OI 6300A line 
and other forbidden lines observed by Hartigan et al., which covered a large sample of T-Tauri 
stars. 

We find that the OI 6300A and the Nell 12.8/im lines are predominantly produced in the 
X-ray-driven wind and thus appear blue-shifted by a few km/s for some of the systems when 
observed at non-edge-on inclinations. We note however that blue-shifts are only produced 
under certain conditions: X-ray luminosity, spectral shape and inner hole size all affect the 
location of the emitting region and the physical conditions in the wind. We caution therefore 
that while a blueshifted line is a tell-tale sign of an outflow, the lack of a blueshift should not be 
necessarily interpreted as a lack of outflow. Comparison of spectrally resolved observations of 
multiple emission lines with detailed model sets, like the one presented here, should provide 
useful diagnostics of the clearing of gaseous discs. 

Key words: accretion, accretion discsxircumstellar matter- planetary systems:protoplanetary 
discs - stars:pre-main sequence 



1 INTRODUCTION 

The circumstellar environments of newly-born low to intermediate 
mass stars have recently been the object of much interest for their 
potential of hosting planetary systems. These (protoplanetary) discs 
of dust and gas hold the reservoir of material from which planets 
may form and provide the medium through which newly formed 
planets and planetesimals migrate towards or away from their par- 
ent star. For this reason no theory of planet formation can be com- 
plete without an understanding of how protoplanetary discs evolve 
and finally disperse. 

The dispersal of protoplanetary discs has been studied obser- 
vationally, through statistics of disc fractions in young clusters (e.g. 
Haisch et al 2001; Hernandez et al 2008; Mamajek 2009) and on 



an object-by-object basis through observations of individual so- 
called transition discs (e.g. Calvet et al 2002), i.e. discs that ap- 
pear to be optically thin in the dust in their inner regions, which 
is generally interpreted as evidence of partial (inside-out) clearing. 
Whether these discs with inner-holes really represent a stage in the 
evolution/dispersal of protoplanetary discs is currently still uncer- 
tain. Currently, the consensus in the literature is that the dispersal 
of discs around T-Tauri stars is characterised by a rather fast tran- 
sition timescale, evidenced by the fact that very few objects (typ- 
ically 10%) are caught in this transition phase (Strom et al. 1989; 
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Strutskie et al. 1990; Ercolano, Clarke & Robitaille 2009; Kim et 
al 2009; Luhman et al 2010; Muzerolle et al 2010fl 

On the basis of these observational results, theoretical efforts 
have thus focused on identifying a mechanism to provide a two- 
timescale dispersal phenomenon, i.e. a mechanism that would al- 
low protoplanetary discs to be detected as optically thick (in the 
dust) for a few million years and then cause them to disperse very 
rapidly in about a tenth of their total lifetime, i.e. a timescale much 
shorter than what would be expected by viscous evolution alone 
(Hartmann et al. 1998, Clarke et al. 2001). 

The proposed mechanisms include planet formation itself 
(e.g. Armitage & Hansen 1999), grain growth (Dullemond & Do- 
minik 2005), photophoresis (Krauss et al 2007), MRI-driven winds 
(Suzuki & Inutsuka 2009), and photoevaporation (Clarke et al 
2001; Alexander et al 2006, Richling & Yorke 2000, Ercolano et 
al 2008b, 2009a, Gorti & Hollenbach 2009; Owen et al. 2010a). 
Understanding which (or what combination) of the proposed mech- 
anisms dominates dispersal is key to understanding disc evolution 
and thus provides important constraints on planet formation the- 
ories by setting the timescale over which planets may form and 
migrate. 

Photoevaporation by soft X-ray radiation (100 eV < E < 
IkeV) has been shown to drive powerful winds (Ercolano et al. 
2008b, 2009a, Gorti & Hollenbach 2009), able to disperse the discs 
rapidly once the accretion rates have fallen to values lower than 
the photoevaporation rates (~10~ 8 M G /yr for X-ray luminosities of 
2xl0 30 erg/s; Owen et al. 2010a; from now on O10). The mass-loss 
rates are roughly proportional to the impinging X-ray luminosity 
(Ercolano 2008b), meaning that the expected mass-loss rates (and 
hence accretion rates) for solar type stars span over the two-three 
decade range of X-rayluminosities(e.g. Preibisch et al. 2005; Giidel 
et al. 2007).The recent observation of shorter disc lifetimes in low- 
metallicity environments (Yasui et al 2009) has been interpreted as 
further evidence for X-ray photoevaporation acting as the main disc 
dispersal mechanism in those clusters (Ercolano & Clarke 2009). 

While encouraging, the evidence in favour of photoevapora- 
tion produced so far consists of statistics derived from observations 
of disc dispersal based on the evolution of their dust component. 
However, dust represents only a percent or so of the total disc mass 
and the evolution of the dust and gas component of a disc do not 
necessarily go hand in hand (e.g. Alexander & Armitage 2007; 
Takeuchi, Clarke & Lin 2004). In the last few years, observations 
of emission lines from warm atomic and low ionisation gas in the 
inner region of discs are starting to provide some helpful clues of 
the physical conditions of the gaseous component in these planet- 
forming regions (Hartigan et al 1995;Herczeg et al 2007; Pascucci 
et al 2007; Razka et al 2007; Espaillat et al 1995; Najita et al 2009; 
Pascucci & Sterzik 2009; Flaccomio et al 2009; Giidel et al 2010, 
submitted). Line emission from low ionisation and atomic species 
is of particular interest, as it is likely to trace material that is lo- 
cated in the warm upper layers of the disc atmosphere or in the 
flow itself, allowing the response of the gas to irradiation from the 
central star to be probed and understood. In particular Pascucci & 
Sterzik (2009) reported evidence of a blue-shifted Nell fine struc- 
ture line at 12.8yum, the blue-shift is consistent with the line being 
formed in a slow moving (a few km/s) disc wind (Alexander 2008). 



There have been claims of longer transition timescales for M-stars or 
older clusters (Sicilia-Aguilar et al. 2008; Currie et al 2009), these results, 
however, are still object of debate in the literature (Ercolano et al 2009, 
Muzerolle et al. 2009, Luhman et al 2010) 



Earlier spectrally resolved observations of optical forbidden lines 
in T-Tauri stars (Hartigan et al., 1995, from now on HEG95), in- 
cluding OI 6300 A, also revealed blue-shifted low velocity com- 
ponents, interpreted as evidence for a disc wind. The interpretation 
of these atomic and low-ionisation emission line profiles is key to 
understand the kinematic structure of the warm upper disc layers 
and its wind. In this paper we provide the tools to exploit currently 
available and future observations using the framework of the X-ray 
driven photoevaporation model of O10 (further discussed in Sec- 
tion[2}. 

Furthermore, with the Herschel Space Telescope and, later, 
James Webb Space Telescope, committed to studying planet for- 
mation through the characterisation of discs, it is both timely and 
important to provide a solid theoretical framework for the inter- 
pretation of these newobservations. In this context this paper has 
two main aims: (i) To provide an atlas of atomic and low ionisa- 
tion emission line intensities and profiles predicted for the photo- 
evaporating protoplanetary discs of O10 at various stages in their 
evolution; (ii) To compare the predictions of OlO's X-ray driven 
photoevaporation model with currently available observations and 
search for the 'smoking gun' of photoevaporation. We therefore 
list line intensities for the brightest lines and calculate line profiles 
for some of the more promising wind/disc tracers. We show re- 
sults for discs at different stages of clearing, i.e. photoevaporating 
primordial discs where an inner hole is yet to form and systems 
in a more advanced clearing stage with inner holes (in dust and 
gas) of increasing radius. We note that, while several authors have 
presented theoretical calculations of emission line luminosities and 
profiles from T-Tauri discs (e.g. Glassgold et al 2007, Meijerink et 
al. 2008, Ercolano et al 2008b, Gorti & Hollenbach 2008, Hollen- 
bach & Gorti 2009, Schisano, Ercolano & Giidel 2010), this is the 
first time that a fully self-consistent radiation-hydrodynamical and 
photoionisation calculation is attempted, which provides us with 
the three-dimensional distribution of gas velocities and emissivi- 
ties in the heated layers of the disc and its wind. 

A comparison of the calculations presented in this paper with 
available observations suggests that X-ray heated discs with hydro- 
dynamically escaping winds produce emission line intensities and 
profiles that are consistent with those that are currently being ob- 
served. It is hoped that the results presented here may be used for 
the interpretation of future observational data which will be crucial 
to further test the validity of the X-ray photoevaporation scenario 
as a dominant disc clearing mechanism. 

The paper is organised as followed. Section 2 briefly describes 
the radiation-hydrodynamic disc models and lists the input param- 
eters. Sections 3 presents the results, including the emission line 
intensities and the predicted profiles for selected transitions. Sec- 
tions 4 and 5 deal with a comparison of our results with currently 
available data. The main points of the paper are briefly summarised 
in Section 6. 



2 THE RADIATION-HYDRODYNAMICS MODELS 

Energetic radiation from a young stellar objects (YSOs) ionises and 
heats the gas in the surrounding protoplanetary disc. The gas ex- 
pands as a result of the extra heating and some regions of the disc 
atmosphere may become unbound, giving rise to a photoevapora- 
tive wind from the disc surface. A calculation of accurate mass loss 
rates and of the density and geometry of the flow requires the simul- 
taneous solution of the two-dimensional radiative transfer and the 
hydrodynamics problem. Owen et al. (2010a, O10) recently pre- 
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sented the first fully self-consistent radiation hydrodynamics cal- 
culations of an X-ray plus EUV photoevaporated protoplanetary 
discs. Ercolano et al (2009) and O10 found irradiation by soft X- 
rays to be most efficient at producing a photoevaporative disc wind 
with mass loss rates of ~ 10~ 8 M o /yr for X-ray luminosities of 10 30 
erg/sec. When coupled with viscous evolution, these wind rates are 
sufficiently large to produce an inner hole in the disc, when the ac- 
cretion rate drops significantly (a factor ~10) below the wind mass- 
loss rates. The remaining material is then dispersed in the final 15- 
20% of the disc lifetime. In this paper we refer to EUV radiation as 
having energies between 13.6eV and 100 eV, soft X-ray radiation 
between 100 eV and 1 keV, and hard X-ray radiation for energies 
larger than 1 keV. The total 'X-ray' luminosity is defined between 
0.1 keV and 10 keV, which is similar to the bandpass quoted by 
Giiedel et al 2009 (0.3-10 keV). 

Two-dimensional hydrodynamic calculations were performed 
using a parameterisation of the gas temperature based on the results 
of three-dimensional radiative transfer and photoionisation calcu- 
lations. Radiative equilibrium was assumed and verified a posteri- 
ori by ensuring that the thermal timescale is shorter than the flow 
timescale. The hydrodynamic calculations were carried out using a 
modified version of the zeus-2d code (Stone & Norman 1992a,b; 
Stone et al 1992) which includes heating from the X-ray and EUV 
radiation from a central YSO. The radiative transfer and photoion- 
isation calculations were carried out using the 3D mocassin code 
(Ercolano et al. 2003, 2005, 2008a) modified according to Ercolano 
et al. (2008b, 2009a). The atomic database included opacity data 
from Verner et al. (1993) and Verner & Yakovlev (1995), energy 
levels, collision strengths and transition probabilities from Version 
5.2 of the CHIANTI database (Landi et al. 2006, and references 
therein) and the improved hydrogen and helium free-bound contin- 
uous emission data of Ercolano & Storey (2006). 

The models presented here are appropriate for a 0.7 M Q star; 
models for lower/higher mass stars are currently being developed 
and will be presented in a future paper. We vary the X-ray lumi- 
nosity to span the observed range for classical T-Tauri stars (cTTs) 
of this mass ( 2xl0 28 <L X < 2xl0 30 erg/sec; e.g. Preibisch et al. 
2005; Gudel et al. 2007) and consider discs at different stages of 
dispersal, from discs extending all the way into their dust destruc- 
tion radius ('primordial discs') to discs with cleared inner holes of 
varying sizes (R,„ = 8.3 AU, 14.2 AU and 30.5 AU, models A, 
D and G in O10, respectively). The ionising spectra used were 
calculated by Ercolano et al (2009a), using the plasma code of 
Kashyap & Drake (2000) from an emission measure distribution 
based on that derived for RS CVn type binaries by Sanz-Forcada 
et al. (2002), which peaks at 10 4 K and fits to Chandra spectra 
of T-Tauri stars by Maggio et al. (2007), which peaks at around 
10 7 5 K. For the remainder of this paper we will refer to the mod- 
els by their X-ray luminosity (0.1 keV < h v < 10 keV), we note, 
however, that the irradiating spectrum has a significant EUV com- 
ponent (13.6 eV < hv < 0.1 keV) with ~L E uv = (see Ercolano et 
al 2009a for further discussion). 

The following elemental abundances were adopted, given as 
number densities with respect to hydrogen: He/H = 0.1, C/H = 
1.4xl0~ 4 , N/H = 8.32 x 10~ 5 , O/H = 3.2 x 10~ 4 , Ne/H = 1.2 x 
10" 4 , Mg/H = 1.1 x 10~ 6 , Si/H = 1.7 x 10~ s , S/H = 2.8 x 10~ 5 . 
These are solar abundances (Asplund et al. 2005) depleted accord- 
ing to Savage & Sembach (1996). D'Alessio et al. (2001) consider 
discs composed of a bimodal dust distribution, where atmospheric 
dust follows the standard MRN model and interior dust consists of 
larger grains with a size distribution still described by a power law 
of index -3.5, but with a„„„ = 0.005 jim and a„ M V = 1 mm. The transi- 



tion between atmospheric and interior dust occurs at a height of 0. 1 
times the midplane gas scale height. As we are mainly interested in 
regions well above the transition point, we have chosen for simplic- 
ity to use atmospheric dust everywhere. Following D'Alessio et al. 
(2001), the dust-to-gas mass ratio of graphite is 0.00025 and that of 
silicates is 0.0004. 

We refer the reader to Ercolano et al. (2008b, 2009a) and O10 
for further details about the codes used and the model setup. 

2.1 Differences in photoevaporation models 

All photoevaporation models are based on the idea that a source of 
radiation (internal - e.g. the YSO, or external e.g. a nearby OB star) 
heats up the gas in the disc, depositing enough energy to leave the 
gas unbound. Low-mass YSOs that are far from external sources 
of energetic radiation, e.g. a typical YSO in Taurus, are the only 
source of irradiation for their discs. A typical YSO produces chro- 
mospheric EUV (13.6 eV < E < 100 eV) and X-ray radiation (E > 
100 eV), with luminosities for the latter spanning several orders of 
magnitude around ~ 10 30 erg/sec. The EUV luminosities are poorly 
known as they cannot be directly observed, but they are estimated 
to be of order 10 41 phot/sec (Alexander et al. 2005). The YSO also 
radiates copiously in the FUV, both from the chromosphere and 
due to accretion shocks (which may also contribute some of the 
X-ray emission). Chromospheric FUV luminosities are compara- 
ble to the observed X-ray luminosities and accretion luminosities 
proportional to the accretion rate of the object. 

Theoretical photoevaporation models that account for one or 
a combination of EUV, X-ray and FUV have been constructed by a 
number of authors and it is worthwhile to point out some of the 
basic differences amongst the current models. In terms of basic 
physics, the models can be split into two categories: (i) Diffuse field 
and (ii) direct field dominated. 

The standard primordial disc EUV models of Hollenbach et 
al (1994), Clarke et al. (2001) and Font et al. (2004) belong to the 
first category. In these models the heating of the upper layers of the 
gaseous disc is dominated by recombination photons from within 
the flow itself. The mass loss rates in this case are accurately deter- 
mined by hydrodynamical calculations and more importantly are 
the mass-loss rate profiles (Font et al. 2004). Mass-loss rates are 
found to be of order 10~ 10 M Q /yr with only a square root depen- 
dence on the ionising luminosity. The mass-loss radial profile is 
narrowly peaked at a few AU for solar type stars. 

As a result of the combined effect of EUV photoevaporation 
and viscous draining, most discs would eventually develop an inner 
hole when the accretion rates fall to a level comparable to the wind 
rates (Clarke et al. 2001, Alexander et al. 2006a). At this point the 
models become direct field dominated (i.e. dominated by radiation 
from the central source) as the inner disc is removed and stellar 
EUV photons can directly reach the inner edge of the photoevapo- 
rating disc. 

The FUV photoevaporation model of Gorti & Hollenbach 
(2009) and the X-ray photoevaporation models of Ercolano et al 
(2008b,2009a) and Owen et al (2010a) also fall in the direct field 
dominated category, both during their primordial and inner hole 
phase. For models in this category, the heating is more or less set 
by the direct flux from the central source. 

Determination of accurate mass-loss rates relies on the abil- 
ity to solve the radiative transfer and thermal calculation problem 
within a hydrodynamical calculation. This has only been done so 
far for the EUV model (Clarke et al 2001, Font et al 2004, Alexan- 
der et al 2006) and the X-ray model (Owen et al 2010), both in 
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the primordial and inner hole phase. A hydrodynamical calculation 
has not been performed for the FUV model and mass-loss rates in 
this case have been estimated by way of an approximation (Gorti 
& Hollenbach 2009), meaning that the values obtained carry poten- 
tially large uncertainties. 

2.1.1 Temperature determination 

The temperature structure in the upper layers of the gaseous disc is 
crucial in determining the flow launching zone and mass loss rates 
(hence the wind structure). The thermal balance of EUV and X- 
ray irradiated models is set by heating by photoionisation against 
cooling which is dominated by collisionally excited line emission. 
The physics and atomic data involved are reasonably well known. 
Ercolano et al (2008b) describes all the heating and cooling mech- 
anisms for the X-ray model in some detail. It is worth noting 
here that Ercolano et al (2008b, 2009a) and Owen et al (2010a) 
use the mocassin code (Ercolano et al 2003, 2005, 2008a) which 
does not include molecular cooling in its thermal balance. The lack 
of molecular cooling has been criticised by Gorti & Hollenbach 
(2009) as a possible source of uncertainty. However, we note that 
the launching region and the associated sonic surface occurs in a 
region where all molecules have already dissociated. While the X- 
rays can deposit energy in the deep molecular layers, their largest 
thermal impact is by far in the upper photodissociated layers. Fur- 
thermore, even at large column densities, where molecular cooling 
can occur, Figure 4 of Glassgold et al (2004) clearly shows that it 
is quickly overwhelmed by dust-gas collisions, which are fully im- 
plemented in mocassin. We therefore conclude that the omission of 
molecular cooling is not a significant source of uncertainty in our 
models (see also further tests and discussion in Owen et al. 2010a). 

In the FUV case explored by Gorti & Hollenbach (2009) the 
heating of the gas is dominated by photoelectric emission from dust 
grains, with PAHs completely dominating the heating channels (see 
e.g. Figure 4 of Gorti & Hollenbach 2008). There are two prob- 
lems with this approach. The first is that the photoelectric yields 
are very poorly known, with uncertainties that could be as large 
as two orders of magnitude (e.g. compare Abbas et al 2006, with 
Feuerbacher & Fitton 1972). Another complication is the fact that 
the abundance (or even presence) of PAH in T-Tauri discs is highly 
uncertain with PAHs detections only in approximately 8% of cases 
(e.g. Geers et al 2006). Indeed models predict that PAHs in the discs 
upper layers are easily destroyed by X-ray irradiation (e.g. Sieben- 
morgen et al 2010). Even assuming that photoelectric yields are 
correct, the uncertainty in the abundance/presence of PAHs in the 
disc makes photoelectric heating almost a free parameter in these 
calculations and the temperature structure largely uncertain. 



3 RESULTS 

The results of our line luminosities and profile calculations (for a 
subset of the lines) are presented in this section. We note that for 
simplicity we choose not to follow the convention of using square 
brackets around forbidden transitions. 

3.1 Atlas of line intensities 

TablesQ]to[4]list the emission line spectra for the models described 
in the previous section. The spectra were calculated using the mo- 
cassin code and using, as the input density structure, the output 
of our modified zeus-2d radiation-hydrodynamic calculations. The 



disk atmosphere is optically thin in the listed lines and therefore 
radiative transfer calculations in the lines is not necessary. Only 
a selected subset of lines is included here, however the full set is 
available upon request from the authors. 

The tables are organised as follows. TablesQ]and[3]list the lu- 
minosities of the collisionally excited transition of heavy elements 
and Tables|2]and|4]list a selection of hydrogen and helium recom- 
bination lines. Tables Q] and [2] show the results for the primordial 
disc models with varying X-ray luminosity (Log(Lx) = 28.3, 29.3 
and 30.3 erg/sec) while Tables[3]and|4]show the results for the inner 
hole models irradiated by Log(Lx) = 30.3 erg/sec and with varying 
hole sizes. 

A casual inspection of these tables immediately reveals that 
the integrated line luminosities are not very sensitive (factors of two 
to five at most) to the size of the inner hole for the range considered 
here (8.3 AU < R,„ < 30.5 AU). However, inner hole models do 
produce, in most cases, collisionally excited lines that are approxi- 
mately a factor two more luminous compared to those of a primor- 
dial disc with the same X-ray luminosity. This is due to the different 
physical conditions in inner hole winds that are less dense and most 
importantly on average warmer than primordial disc winds. 

Line luminosities are, in general, correlated with X-ray lumi- 
nosity, the response of different lines to changes in L x is however, 
different for different transitions, and depends on how sensitive to 
changes in luminosity a line's emitting region is, and on the depen- 
dance of a given emission line to temperature, density or ionisation 
structure. A number of factors are relevant to determine which re- 
gion contributes to the emission of a given line. The ionic abun- 
dance is of course important; collisionally excited lines of Nell and 
OI require Ne + and O , as well as colliding particles, like (e.g.) 
electrons and/or neutral hydrogen. HI recombination lines arise in 
region where there is a non zero abundance of H + and e - . Another 
key property is temperature. Collisionally excited lines depend ex- 
ponentially on the gas temperature through the Boltzmann term and 
therefore preferentially trace the hotter regions. In the case of col- 
lisionally excited lines in the infrared, however, this temperature 
dependence is weakened by the fact that the gas temperature is gen- 
erally much higher than the line excitation temperature, causing the 
Boltzmann term to tend to unity. Recombination lines on the other 
hand, only have an inverse square root dependence on temperature 
and therefore are emitted preferentially in cooler regions, as long 
as these are not completely neutral. Gas density also influences the 
emissivity of collisionally excited lines through critical density ar- 
guments. 

It should also be noted that increasing the X-ray luminosity 
not only affects the ionisation rate of the wind, but also affects its 
density. Higher luminosities drive a more vigorous (denser) wind, 
with a roughly linear correlation between L* and gas density in 
the wind. For the primordial discs presented here the scaling of 
the density with L x can be understood simply. Given the wind is 
thermally driven, the specific energy requirements to drive the flow 
will only depend on the local value of the effective potential. Once 
a freely expanding wind is established it will remain at approx- 
imately constant temperature. Since for transonic/sonic flows the 
velocity is strongly coupled to the sound speed of the gas, which 
in turn only depends on the value of the gas temperature, then, for 
all density and temperature regimes we are interested in, this lo- 
cal specific energy requirement is equivalent to a local temperature 
condition in the launch regiorQ. Thus the spatial variation of tem- 



2 We have and will implicitly assume that the total input energy from the 
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Figure 1. Line profiles for (from the top) Nell 12.8^m, OI 6300A and a generic hydrogen recombination line. The left, middle and right panels are for 

/■40 i 

primordial disc models irradiated by Log(Lx) = 28.3, 29.3 and 30.3, respectively. These line profiles are all normalised, such that J 4() L[i'/(km s )] dv = 1. 
The profiles are colour-coded according to their inclination angle, a colour legend in given in the top-left panel. 



perature within the flow is approximately set by the mass of the 
central star, and will be independent of the local density and X- 
ray luminosity. Given also that temperature can be expressed as a 
monotonic function of the ionisation parameter, the spatial varia- 
tion of the ionisation parameter is approximately fixed by the mass 
of the central star. Hence, using the definition of the ionisation pa- 
rameter £ =L x /nr 2 and given the value of the ionisation parameter 
is spatially fixed, at a point r =constan£| we find that gas density 
n ocLjo as can be seen in the primordial discs shown in Figure[3] 

An element-by-element discussion, comments on observabil- 
ity of the various transitions and a comparison with observations 
and other theoretical investigation is postponed to Section 4 



high energy radiation is sufficient to balance the work done driving any flow. 
See O10 of a detailed discussion of this condition in relation to X-ray driven 
photoevaporative winds. 

3 Note throughout this paper we use (r, 6, <f>) and (R, ip, z) to distinguish be- 
tween spherical and cylindrical polar co-ordinates respectfully. 



3.2 Line profiles 

We have used our two-dimensional map of emissivities and gas 
velocities to reconstruct a three-dimensional cube of the disc and 
calculate the line-of-sight emission profiles for a number of tran- 
sitions. The emissivities are calculated using 120 radial points 
(N R = 120) and 1500 height points (N- = 1500). By assum- 
ing azimuthal symmetry about the Z-axis, and reflection symme- 
try about the disc mid-plane the resultant 3D grid has dimensions 
N R x 2N- x A^. We adopted a value of N v = 130 and checked that 
this was large enough to fully resolve the profiles. For the chosen 
emission lines the disc atmosphere is optically thin and therefore 
we assume that the line contribution of each parcel of gas can es- 
cape freely along the line of sight, provided that this does not in- 
tercept the midplane (e.g. Alexander 2008, Schisano et al. 2010). 
For almost edge-on inclinations escape columns are large enough 
for dust attenuation to start playing a role. This is not accounted 
for in this work and the magnitude of this effect depends largely on 
poorly known parameters such as dust to gas ratios, the degree of 
dust settling and grain growth in the disc. 
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Figure 2. Line profiles for (from the top) Nell 12.8/jm, OI 6300A and a generic hydrogen recombination line. The left, middle and right panels are for inner- 
hole disc models irradiated by Log(Lx) = 30.3, with inner holes of radius 8.3, 14.2 and 30.5, respectively. These line profiles are all normalised, such that 
f40 _i 

J 40 L[\'/(km s )] dv = 1. The profiles are colour-coded according to their inclination angle, a colour legend in given in the top-left panel. 



The line luminosity is then computed by including a Doppler 
broadening term in each cell. Hence the line luminosity at a given 
velocity v is computed by evaluating the following integral by di- 
rect summation. 



L(V): 



y/2nv,i,(r) : 



: exp - 



[ v - Vi„s(r)Y 
2v th {rY 



(1) 



Where E{r) is the volume averaged power emitted at a point r, v; os (r) 
is the projected gas velocity along the line of sight and v th is the 
local thermal velocity (rms velocity) of the emitting atom given by 



3K B T(r) 

fliin H 



(2) 



With Hi the mean molecular weight of the required atom. The reso- 
lution of the velocity array is 0.25 km/s. 

We have evaluated the profiles of the OI 6300A, the 
Nell 12.8/um and a generic HI recombination line, as these have 
been the focus of recent observational studies (e.g. HEG95, Pas- 
cucci & Sterzik 2009; Najita et al 2009). Profiles for any of the 



lines listed in the table (and many others that are not) can easily be 
calculated upon request by the authors. 

Figure Q] shows the normalised profiles of the Nell fine struc- 
ture line at 12.8yum, the OI forbidden line at 6300 A and a hy- 
drogen recombination line for the primordial disc model with 
Log(L x ) = 28.3 erg/s (left), 29.3 erg/s (middle) and 30.3 erg/s 
(right). Figure [2] shows the same line profiles for inner-hole mod- 
els with Log(L^) = 30.3 erg/s and inner hole radii of 8.3 AU (left), 
14.2 AU (middle) and 30.5 AU (right). Profiles are computed for 
ten inclination angles from 0" to 90° degrees in steps of 10°. The 
profiles were degraded to an instrumental resolution of R = 100000, 
typical of current ground based facilities. 

Figure[3]shows 2D density maps (colour gradients) of the pri- 
mordial disc (left panels) and inner-hole models (right panels), re- 
spectively, with velocity vectors overlain. The models from top to 
bottom are for increasing X-ray luminosities (left) and increasing 
inner hole radii (right). The black, magenta and yellow contours 
delineate the region containing 85% of the emission for the HI re- 
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Figure 3. Density maps showing the location of the 85% emission region of HI recombination lines (black contour), OI 6300A line (yellow contour) and 
Nell 12.8 fim line (magenta contour). The velocity field is represented by the red arrows. Left: Primordial disc models irradiated by Log(Lx) = 28.3 erg/sec 
(top), Log(Lx) = 29.3 erg/sec (middle) and Log(Lj-) = 30.3 erg/sec (bottom). Right: Inner hole models irradiated by Log(L^) = 30.3 erg/sec with inner hole 
radii of R,„ = 8.3 AU, R,„ = 14.2 AU and R,„ = 30.5 AU 
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Table 1. Collisionally excited lines. Primordial disk models with varying X-ray luminosity. Wavelengths are in vacuum. Only selected lines are listed here, the 
full set is available upon request from the authors. 



Species 


Wavelength 


Line Luminosity [L Q 


1 


Species 


Wavelength 


Line Luminosity [L Q 






[A] 


Lo2("Lv)=28 3 


29.3 


30.3 




[A] 


LosfLvi=28 3 


29.3 


30.3 


CI 


8729.51 


1.79E-08 


8.21E-08 


5.78E-07 


Mgl 


4573.180 


6.11E-08 


8.21E-07 


1.08E-05 


CI 


9826.85 


4.32E-08 


3.67E-08 


6.23E-07 


Mgl 


4564.670 


1.37E-09 


4.73E-09 


2.16E-08 


CI 


9852.99 


1.27E-07 


1.06E-08 


1.16E-06 


Mgl 


2852.140 


8.05E-09 


1.16E-07 


1.54E-06 


CII 


1577287. 


8.80E-08 


2.80E-07 


6.64E-07 


Mgll 


2803.530 


4.77E-09 


7.92E-08 


1.28E-06 


CII 


2325.400 


3.13E-09 


5.60E-08 


9.43E-07 


Mgll 


2796.350 


9.40E-09 


1.54E-07 


2.47E-06 


CII 


2324.210 


3.13E-09 


5.02E-08 


5.87E-07 


Sin 


348189.4 


5.67E-07 


1.06E-06 


1.85E-06 


CII 


1334.530 


2.08E-10 


4.35E-09 


7.42E-08 


Sill 


2335.120 


1.20E-10 


2.03E-09 


3.24E-08 


CII 


2328.840 


3.73E-09 


6.67E-08 


1.11E-06 


Sill 


1808.010 


1.43E-10 


2.61E-09 


4.17E-08 


CII 


2327.640 


9.84E-09 


1.54E-07 


1.85E-06 


Sill 


2350.890 


7.90E-11 


1.35E-09 


2.00E-08 


CII 


2326.110 


2.08E-08 


3.57E-07 


5.56E-06 


Sin 


2344.920 


3.87E-10 


6.57E-09 


1 .00E-07 


CII 


1335.660 


4.02E-11 


8.60E-10 


1.48E-08 


sin 


2335.320 


5.52E-10 


9.57E-09 


1 .46E-07 


CII 


1335.710 


3.58E-10 


7.83E-09 


1.34E-07 


Sill 


1816.930 


2.53E-10 


4.64E-09 


7.26E-08 


OI 


6302.030 


2.83E-07 


1.93E-06 


1.25E-05 


SII 


6732.690 


2.08E-07 


1.16E-06 


1.08E-05 


OI 


6365.530 


9.10E-08 


6.28E-07 


4.01E-06 


SII 


6718.310 


1.41E-07 


7.44E-07 


5.72E-06 


OI 


2973.160 


6.11E-10 


9.18E-09 


1.06E-07 


SII 


4077.510 


5.81E-08 


4.25E-07 


3.24E-06 


OI 


5578.890 


5.22E-09 


7.83E-08 


9.12E-07 


SII 


4069.760 


2.23E-07 


1.64E-06 


1.22E-05 


Oil 


3729.880 


1.79E-08 


1.06E-07 


6.49E-07 


SII 


10339.23 


3.87E-08 


2.80E-07 


2.16E-06 


Oil 


3727.090 


3.58E-08 


2.51E-07 


1.85E-06 


SII 


10289.55 


5.52E-08 


4.06E-07 


2.93E-06 


on 


2470.970 


2.23E-09 


2.70E-08 


2.78E-07 


SII 


10373.34 


1.64E-08 


1.16E-07 


9.12E-07 


on 


2471.090 


6.86E-09 


8.70E-08 


8.65E-07 


SII 


10323.33 


6.56E-08 


4.83E-07 


3.55E-06 


on 


7320.910 


1 . /9E-U9 


1 i OT7 no 

1. 12E-U8 


2.16E-07 


SIII 


335008.4 


1.64L-08 


1.25E-0 7 


1.28E-06 


Oil 


7321.980 


4.17E-09 


5.41E-08 


5.41E-07 


sin 


187055.7 


2.38E-08 


1.64E-07 


2.00E-06 


Oil 


7331.700 


2.83E-09 


3.57E-08 


3.55E-07 


SIII 


9070.050 


3.58E-09 


5.41E-08 


8.19E-07 


on 


7332.780 


2.23E-09 


2.90E-08 


2.93E-07 


SIII 


3722.450 


2.23E-10 


5.02E-09 


8.34E-08 


Nell 


128155.8 


9.25E-08 


5.51E-07 


5.41E-06 


SIII 


9532.250 


2.08E-08 


3.09E-07 


4.63E-06 


Nelll 


155545.2 


1.01E-08 


1.16E-07 


1.06E-06 


SIII 


1728.950 


2.23E-11 


7.73E-10 


2.16E-08 


Nelll 


3869.850 


2.98E-09 


6.28E-08 


8.34E-07 


SIII 


6313.650 


4.02E-10 


8.89E-09 


1.48E-07 


Nelll 


360230.6 


5.96E-10 


6.47E-09 


5.72E-08 


SIV 


105108.3 


4.92E-08 


2.99E-07 


3.71E-06 


Nelll 


3968.580 


8.80E-10 


1.93E-08 


2.47E-07 














-40 -20 20 -20 20 -20 20 -20 20 40 

Velocity [km s" 1 ] 

Figure 4. Representative shapes for the line profile shapes calculated from our models. Profile A (left) can be well approximated by a Gaussian. Profile B is 
a Gaussian with large wings. Profile C is a double peaked profile and profile D are profiles where two clear peaks are not resolved, yet the profile is broad, 
flat-topped or significantly asymmetric. The blue vertical dots show the location of v„ eB t as given in Table [5] and the red horizontal dashed line shows the 
FWHM for the given profiles. 



combination lines, the Nell 12.8yurr0 and the OI 6300A, respec- 
tively. The noisy appearance of some of the emitting regions are 
due to Monte Carlo fluctuations and the imposition of a hard cut- 
off to draw the contours. We have checked that the noise does not 
affect luminosities and line profiles to more than a 5% level. 

The OI 6300A line (yellow contours) originates almost exclu- 
sively in the wind or in the very inner (hot and bound-r? < 1 AU) 

4 The Nelll 15.5/im region is very similar. 



region of the disc (in the primordial case), which is where the tem- 
peratures are high enough to populate its upper level. The same 
would be true for all collisionally excited optical and UV lines. It 
is the contribution from the bound hot inner disc that produces the 
'winged' profiles for the OI 6300A clearly seen in the lower lumi- 
nosity primordial models. 

The Nell mid-infrared (MIR) lines can arise in much cooler 
regions and therefore, while still predominantly produced in the 
wind, they are not confined to the hot inner layers, but span over a 
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Table 3. Collisionally excited lines. Inner hole disk models with varying hole size. Wavelengths are in vacuum. Only selected lines are listed here, the full set 
is available upon request from the authors. 



Species 


Wavelength 


Line Luminosity [L, 


3 ] 


Species 


Wavelength 


Line Luminosity [L, 


•>] 




[A] 


R,„=8.3 AU 


14.2 AU 


30.5 AU 




[A] 


R,„=8.3 AU 


14.2 AU 


30.5 AU 


CI 


8729.51 


5.57E-07 


5.57E-07 


5.38E-07 


Nelll 


3869.850 


2.25E-06 


1.99E-06 


2.61E-06 


CI 


9826.85 


7.52E-07 


5.57E-07 


7.02E-07 


Nelll 


360230.6 


1.21E-07 


1.57E-07 


2.24E-07 


CI 


9852.99 


1.30E-06 


1.66E-06 


2.09E-06 


Nelll 


3968.580 


6.86E-07 


6.14E-07 


7.48E-07 


CII 


1577287. 


1.12E-06 


1.41E-06 


2.09E-06 


NelV 


2422.360 


1.56E-08 


1.90E-09 


8.98E-08 


CII 


2325.400 


5.03E-07 


4.48E-07 


3.96E-07 


Mgl 


4573.180 


1.12E-06 


9.13E-07 


6.66E-07 


CII 


2324.210 


5.38E-07 


4.81E-07 


4.26E-07 


Mgl 


4564.670 


5.99E-08 


7.22E-08 


9.72E-08 


CII 


1334.530 


1.04E-07 


9.13E-08 


8.98E-08 


Mgl 


2852.140 


2.69E-07 


2.32E-07 


1.87E-07 


CII 


2328.840 


5.99E-07 


5.31E-07 


4.71E-07 


Mgll 


2803.530 


1.04E-06 


9.13E-07 


7.48E-07 


CII 


2327.640 


1.65E-06 


1.49E-06 


1.27E-06 


Mgll 


2796.350 


1.99E-06 


1.82E-06 


1.49E-06 


CII 


2326.110 


3.30E-06 


2.90E-06 


2.61E-06 


Sill 


348189.4 


1.65E-06 


1.66E-06 


1.87E-06 


CII 


1335.660 


1.99E-08 


1.74E-08 


1.79E-08 


Sin 


2335.120 


1.99E-08 


1.82E-08 


1.64E-08 


CII 


1335.710 


1.82E-07 


1.57E-07 


1.64E-07 


Sill 


1808.010 


4.08E-08 


3.65E-08 


3.59E-08 


OI 


6302.030 


4.51E-05 


4.48E-05 


3.59E-05 


Sill 


2350.890 


1.30E-08 


1.16E-08 


1.12E-08 


OI 


6365.530 


1.39E-05 


1.41E-05 


1.12E-05 


Sill 


2344.920 


6.42E-08 


5.72E-08 


5.31E-08 


OI 


2973.160 


1.73E-08 


1.07E-08 


5.98E-09 


Sill 


2335.320 


9.55E-08 


8.30E-08 


7.48E-08 


OI 


5578.890 


1.47E-07 


9.13E-08 


5.16E-08 


Sill 


1816.930 


7.21E-08 


6.39E-08 


6.28E-08 


on 


3729.880 


1.91E-06 


2.32E-06 


3.44E-06 


SII 


6732.690 


2.78E-05 


3.07E-05 


3.14E-05 


on 


3727.090 


5.47E-06 


6.64E-06 


8.98E-06 


SII 


6718.310 


1.47E-05 


1.57E-05 


1.87E-05 


on 


2470.970 


9.55E-07 


9.13E-07 


7.40E-07 


SII 


4077.510 


7.73E-06 


6.64E-06 


5.01E-06 


on 


2471.090 


3.30E-06 


3.15E-06 


2.76E-06 


SII 


4069.760 


2.43E-05 


1.99E-05 


1.49E-05 


on 


7320.910 


7.73E-07 


6.97E-07 


5.91E-07 


SII 


10339.23 


5.12E-06 


4.40E-06 


3.29E-06 


on 


7321.980 


2.08E-06 


1.99E-06 


1.72E-06 


SII 


10289.55 


5.90E-06 


4.89E-06 


3.59E-06 


on 


7331.700 


1.30E-06 


1.16E-06 


9.72E-07 


SII 


10373.34 


2.08E-06 


1.82E-06 


1.34E-06 


on 


7332.780 


1.12E-06 


1.07E-06 


8.98E-07 


SII 


10323.33 


6.86E-06 


5.72E-06 


4.19E-06 


OIII 


883392.2 


/.U3L-U8 


/. /2L-08 


1.12E-07 


sin 


335008.4 


t.8zr,-Uo 


1 OOT7 (\C 


2.32E-06 


OIII 


518134.7 


2.78E-07 


3.23E-07 


4.71E-07 


sin 


8830.960 


1.65E-08 


1.41E-08 


1.49E-08 


OIII 


4960.290 


5.38E-07 


3.98E-07 


1.87E-06 


SHI 


187055.7 


3.38E-06 


3.65E-06 


4.04E-06 


OIII 


2321.670 


1.39E-08 


6.22E-09 


6.88E-08 


sin 


9070.050 


1.91E-06 


1.66E-06 


1.79E-06 


OIII 


5008.240 


1.65E-06 


1.16E-06 


5.53E-06 


sin 


3722.450 


1.30E-07 


1.07E-07 


1.49E-07 


OIII 


1666.150 


1.91E-08 


7.30E-09 


1.12E-07 


sin 


9532.250 


1.12E-05 


9.13E-06 


1.04E-05 


OIII 


4364.450 


1.56E-08 


7.05E-09 


7.48E-08 


sin 


1728.950 


2.51E-08 


1.99E-08 


4.86E-08 


Nell 


128155.8 


1.21E-05 


1.16E-05 


1.04E-05 


sin 


6313.650 


2.34E-07 


1.82E-07 


2.61E-07 


Nelll 


155545.2 


2.43E-06 


2.90E-06 


3.44E-06 


SIV 


105108.3 


5.73E-06 


5.64E-06 


5.08E-06 



larger region. The Nell emitting region (magenta contours in Fig- 
ure [3} and line profiles vary significantly with X-ray luminosity, 
tracing the highest electron density regions (L(NeII) ocnj). 

The profiles of HI recombination lines are more complicated 
as they have contributions arising both from the wind and from 
the bound disc. However, as will be discussed later, our models 
underpredict the luminosities of observed IR hydrogen recombina- 
tion lines by one or two orders of magnitude (Pascucci et al 2007; 
Ratzka et al 2007; Najita et al 2010) and underpredict the lumi- 
nosities of observed Pa/3 and Bry lines by four orders of magnitude 
(Muzerolle et al, 1998), implying that the emission must be dom- 
inated by regions that are not included in our models (e.g. dense 
plasma close to the star). The resulting line profiles may therefore 
bear little resemblance to those shown in this paper. The left panels 
of Figure[3](black contours) show that in the primordial disc models 
there is only a small wind contribution, which results in the hydro- 
gen recombination profiles showing no blue-shift atall. While inner 
hole models have larger wind contributions to the HI recombination 
lines, the resulting profiles are still not blue-shifted, although this 
is virtually true for all line profiles from inner hole models, even 
those formed exclusively in the wind. 

The reason why inner hole models produce very small or no 
blue-shifts is that at low inclinations both the blue and red-shifted 
components of the wind are visible through the hole, resulting in 



only a small, often undetectable, blue-shift. At higher inclination 
the lines are Keplerian and no blue-shift is expected anyway. We 
will discuss the significance of this result in the next section, when 
comparisons are made with observations, but we note at this point 
that this is only true if the winds are optically thin to the lines in 
question. 



To aid comparison with observations we have summarised the 
main profile characteristics in Tables [5] and [6] These include the 
velocity of the peak(s), v pe „ k , the full-width-half-maxima (FWHM) 
and the profile "type". We have classified the profiles with four 
main types as shown in Figure [4] Type "A" is a single peaked pro- 
file that could be well approximated by a Gaussian. Type "B" is a 
single peaked Gaussian-like profile, but with large wings. This is 
generally produced when there is both a wind and bound disc con- 
tribution to the line and/or the emission region is spread over a large 
radial domain. Type "C" is a double peaked profile, in this case the 
two values of v pm k are given. Type "D" are profiles where two clear 
peaks are not resolved, yet the profile is broad, flat-topped or sig- 
nificantly asymmetric; in this case the value at the centre of the line 
(i.e. the middle of the FWHM segment, horizontal dashed red line 
in Figure[4]l is given; strictly speaking, this is not the location of the 
'peak' however we keep the label v pea k for simplicity. 
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Table 2. Recombination Lines of H, Hel and Hell. Primordial disc models 
with varying X-ray luminosity. Wavelengths are in vacuum. Only selected 
lines are listed here, the full set is available upon request from the authors. 



Hydrogen 



Xr3.n s iti o n 


Wavelength 


Line Luminosity [L Q 


1 






Lo°(Ly)=28 3 


29.3 


30.3 






HI 






3-2 


6564.696 


2.62E-07 


2.77E-06 


1.10E-05 


4-2 


4862.738 


7.87E-08 


9.18E-07 


3.20E-06 


5-2 


4341.730 


3.52E-08 


4.26E-07 


1.43E-06 


6-2 


4102.935 


1.94E-08 


2.41E-07 


7.97E-07 


5-3 


12821.67 


1.62E-08 


1.67E-07 


6.70E-07 


6-3 


10941.16 


8.38E-09 


9.13E-08 


3.40E-07 


7-3 


10052.19 


5.03E-09 


5.70E-08 


2.04E-07 


8-3 


9548.649 


3.35E-09 


3.91E-08 


1.35E-07 


7-4 


21661.29 


3.05E-09 


3.01E-08 


1.26E-07 


7-5 


46537.93 


2.20E-09 


1.88E-08 


9.64E-08 


7-6 


123719.3 


1.68E-09 


1.23E-08 


7.99E-08 


8-7 


190619.3 


7.79E-10 


5.50E-09 


3.76E-08 


9-7 


113087.2 


5.07E-10 


3.86E-09 


2.35E-08 






Hel 






4 3 D-2 3 P 


4473.000 


8.17E-10 


1.63E-08 


3.50E-08 


3 3 P-2 3 S 


3890.000 


1.81E-09 


3.79E-08 


8.50E-08 


3 3 D-2 3 P 


5877.000 


2.38E-09 


4.79E-08 


1.03E-07 


2 3 P-2 3 S 


10833.00 


1.82E-08 


4.43E-07 


1.10E-06 






Hell 






4-3 


4689.069 


1.77e-10 


3.91e-09 


5.92e-08 



Table 4. Recombination lines of H, Hel and Hell. Inner hole disc models 
with varying hole size. Wavelengths are in vacuum. Only selected lines are 
listed here, the full set is available upon request from the authors. 



Transition Wavelength Line Luminosity [L Q ] 





[A] 


8.3 


14.2 


30.5 


3-2 


6564.696 


1.40E-05 


2.43E-05 


4.76E-06 


4-2 


4862.738 


4.58E-06 


7.88E-06 


1.55E-06 


5-2 


4341.730 


2.10E-06 


3.60E-06 


7.11E-07 


6-2 


4102.935 


1.15E-06 


1.97E-06 


3.90E-07 


5-3 


12821.67 


8.28E-07 


1.44E-06 


2.81E-07 


6-3 


10941.16 


4.42E-07 


7.67E-07 


1.49E-07 


7-3 


10052.19 


2.67E-07 


4.61E-07 


9.03E-08 


8-3 


9548.649 


1.75E-07 


3.02E-07 


5.90E-08 


7-4 


21661.29 


1.47E-07 


2.59E-07 


5.02E-08 


7-5 


46537.93 


9.62E-08 


1.71E-07 


3.31E-08 


7-6 


123719.3 


6.62E-08 


1.20E-07 


2.30E-08 


8-7 


190619.3 


2.98E-08 


5.42E-08 


1.03E-O8 


9-7 


113087.2 


2.04E-08 


3.70E-08 


7.12E-09 






Hel 






4 3 D-2 3 P 
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3.3 Nell 12.8yum and OI 6300A radial profiles 

The model's radial intensity profiles of the Nell 12.8/um and of the 
OI 6300A line are shown in the left and right panels of Figure [331 
The top panels show the response of the radial intensity profiles 
of these lines to changes in the X-ray luminosity. Wider profiles 
are expected for lower X-ray luminosities since these drive a wind 



that is less dense than in the higher X-ray luminosity case, mean- 
ing that the emission region is more extended. The bottom panels 
of Figure [331 show the radial intensity profiles for discs irradiated 
by the same X-ray luminosity (2xl0 30 erg/sec), with the black line 
representing a primordial disc and the red, green and blue lines rep- 
resenting discs with inner holes of size 8.3, 14.2 and 30.5 AU, re- 
spectively. Both the Nell and OI lines show a very clear sequence 
with inner hole size suggesting that spectro-astrometry (e.g. Pon- 
toppidan et al 2009) may perhaps be used to test the evolutionary 
states of gaseous discs in the clearing phase. We draw attention to 
the fact that the peak of the emission in all inner-hole cases is at 
radii smaller than the dust inner radius of the disc. This is obvious 
by looking at the structure of the wind and the location of the emis- 
sion region (see Figure[3j- So the presence of gas in the inner hole 
of these objects is simply due to gas that is being photoevaporated 
from the inner edge of the disc. 

3.4 Caveats and limitations 

Before moving on to a comparison of our results with observations 
it is worth stating a few caveats and limitations of our models. 

• It is important to realise that the current models can only pro- 
vide a snapshot of the various flow solutions and line profiles that 
will arise. This means that it is currently not possible to present any 
form of 'interpolation' scheme for the line profiles and only general 
trends can be qualitatively extracted. There are two main reasons 
for this: (i) Most of the emitting regions (see Figure [3} is located 
in the subsonic region of the flow. In this region the flow veloc- 
ity and direction is strongly and non-linearly affected by the inter- 
action between the centrifugal force term and velocity divergence 
along the stream line (Ercolano et al. 2009a). Hence the variation 
of the blue shift with inclination and X-ray luminosity can only 
be understood in a qualitative way. The lack of a simple relation 
between line luminosity and inclination or X-ray luminosity is par- 
ticularly evident in the behaviour of the OI line, (ii) The emitting 
region that dominates the luminosity for a given line can be con- 
sidered a function of three variable: temperature, gas density and 
electron/ion fraction. All three influence the emissivity of a given 
pocket of gas and the effect they have on each other again arises 
from highly non-linear phenomena, meaning the emission region 
can change smoothly with X-ray luminosity (e.g. HI) or smoothly 
for some values of X-ray luminosity (e.g. h x > Mr erg/s) and 
rapidly for others (e.g. h x < 10 29 erg/s), as in the Nell line. 

• The ionising spectrum, particularly in the soft X-ray and EUV 
range, is somewhat uncertain. Observations of this part of the spec- 
trum are hindered by absorption from interstellar hydrogen along 
the line of sight or closer to the star. While we have attempted to 
use a realistic spectrum (see Section[2]and discussion in Ercolano et 
al. 2009a), it is very likely that T-Tauri stars posses a range of spec- 
tra of different hardness. Accretion funnels close to the star may 
also screen the disc from EUV and soft X-ray effectively changing 
the hardness of the radiation field seen by the disc. The effects of 
screening on photoevaporation rates have been explored and dis- 
cussed in detail by Ercolano et al (2009a). 

• Our models assume azimuthal symmetry, however the effects 
of (possibly inclined) magnetic fields on corotation radius may 
cause twists and warps at the disc's inner edge (e.g. Bouvier et al 
2007; Romanova et al 2003; Long et al 2007). Furthermore accre- 
tion funnels may only screen a fraction of the radiation resulting in 
azimuthally varying irradiation. As the structure of magnetic fields 
around T-Tauri stars is becoming clearer through observations and 
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Radius [AU] Radius [AU] 

Figure 5. Radial intensity profiles of Nell 12.8/mi (left panels) and OI 6300A (right panel). The top panels show how the profiles respond to variations in the 
X-ray luminosity in the case of a primordial disc, while the bottom panels show the response to different hole sizes (R,„). R su b is the dust sublimation (destruction) 
radius. All models shown in the bottom panels are irradiated by an X-ray luminosity of 2xl0 30 erg/sec. 



modeling (Donati et al 2007, 2008; Gregory et al 2008; Jardine 
2008; Hussain et al 2009; Skelly et al 2009), it will be very inter- 
esting to test these effects on photoevaporation and disc emission. 

• The models presented in this paper are for a 0.7 M Q star. It 
is difficult to extrapolate these results to different mass stars, how- 
ever models appropriate for different mass stars are currently being 
developed (Owen et al 2010b, in prep). 

• Our models do not account for molecular chemistry. While 
this is unlikely to produce a significant error in the photoevapora- 
tion rates and wind structure (see discussion in O10), it may affect 
the luminosities of low excitation temperature emission lines from 
neutral species (e.g. fine structure lines of OI and CI), which are 
produced in a region close to the atomic/molecular transition (the 
photodissociation region, PDR). We note that all the lines listed in 
Tables [TJ to [4] however, are not significantly affected by the PDR, 
due to their high excitation temperatures (optical/UV collisionally 
excited lines) and/or ionisation requirements. 

• Excitation rates of Ne + by collisions with neutrals are not 
known, therefore theoretical calculations of this line, including this 



work, only account for collisions with electrons. This is likely to 
result in the Nell 12.8/jm line being underestimated by a factor of 
~2 (Glassgold et al 2007). 

• The underlying disc flaring could be important since it deter- 
mines the solid angle subtended by the disc to the ionising radia- 
tion. Owen et al (2010) briefly discuss this point with regards to 
mass loss rates in their Section 5, point vii. A discussion of the 
effects of flaring on line luminosities, specifically for the case of 
the Nell 12.8/jm and Nelll 15.3yum lines, is given by Schisano, Er- 
colano & Giidel (2010). 



4 DISCUSSION: THE EMISSION LINE SPECTRA 

The general behaviour of emission line intensities with X-ray lu- 
minosity and inner hole size has already been discussion in Sec- 
tion |3.1| Here we further comment on an element-by-element basis 
on the significance and observability of some of the transitions. 
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Table 5. Line profiles from primordial discs irradiated by Log(Lx) = 28.3, 29.3 and 30.3 erg/s. The table lists the velocity of the peak, the full-width-half- 
maximum (FWHM) and the type of the profile (see Figure 4 for details), for ten inclinations from 0° to 90°. The profiles were degraded to an instrumental 
resolution of R = 100000. 
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4.1 Carbon 

CI fine-structure lines are produced in a region that is very close 
to the atomic to molecular transition. The luminosities calculated 
by mocassin for these transitions are likely to carry some error and 
have therefore been omitted from the table and will be discussed 
in a forthcoming paper where a simplified chemical network is in- 
cluded. These lines may be interesting as probes of the gas at larger 
radii (e.g. Ercolano et al 2009b). 

The CII fine structure line at 157/zm, however, should be less 
affected by the molecular/atomic transition and therefore we dis- 
cuss it further here. We report a luminosity of ~ 10~ 6 L Q for an X-ray 
luminosity of 2xl0 30 erg/sec, this should make this line an interest- 
ing target for Herschel observations. We note that the luminosity we 
predict for this line is almost three orders of magnitude larger than 



what was calculated by Meijerink et al. (2008); as discussed in Er- 
colano et al (2008b), the low critical density for this line (49 cm -3 ) 
implies that it will be preferentially emitted in low density regions, 
such as the upper disc atmosphere and the flow, which were not 
included in the Meijerink calculations. 

We also list a number of UV lines from CII, which are emitted 
predominantly in the wind; the strongest of these transitions is the 
doublet at 2327 A. 



4.2 Oxygen 

As in the case of carbon, we limit our discussion of neutral oxygen 
emission to forbidden lines with high temperature Boltzmann terms 
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Table 6. Line profiles from inner-hole discs irradiated by Log(Lx = 30.3 erg/s, with inner hole radii of 8.3, 14.2 and 30.5 AU. The table lists the velocity of 
the peak, the full-width-half-maximum (FWHM) and the type of the profile (see Figure 4 for details), for ten inclinations from 0° to 90°. The profiles were 
degraded to an instrumental resolution of R = 100000. 
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and exclude the 01 fine-structure lines that may be uncertain in our 
calculations for the reasons discussed above. 

The 01 6300A forbidden line deserves special attention. This 
line was detected by HEG95 in the spectra of several T-Tauri stars. 
Detection of this line in T-Tauri stars has been reported by vari- 
ous other authors including White & Hillenbrand (2004) and Mo- 
hanty et al (2005). The OI 6300A line profiles reported by HEG95 
all showed a double peak, with one component blue-shifted to a 
few hundred km/s (high velocity component HVC) and a second 
component typically blue-shifted by only a few km/s (low velocity 
component, LVC). The HVC is generally attributed to spatially un- 
resolved emission from a dense jet close to the star. The LVC, seen 
also in other forbidden lines, such as OI 5577 A, and SII 6731 A, 
was interpreted by these authors as being produced in a slow mov- 



ing disc wind. This was evidenced by the fact that blue-shifts were 
largest for forbidden lines of lowest critical density, suggesting 
that the low velocity component accelerates away from the surface 
of the disc. HEG95 also noted that all optically thick discs that 
they observed had LVC forbidden line emission detections. Fur- 
thermore, all the LVCs were blue-shifted to within the measure- 
ments errors (±2 km/s), leading them to conclude that, if the LVC 
arises from a disc wind, then all objects with optically thick inner 
discs must have a disc wind. 

While HEG95 did not ascertain the nature of the disc wind, a 
photoevaporative wind origin was investigated by Font et al (2004), 
who explored the ability of an EUV-driven photoevaporative wind 
to explain the LVC lines observed in the HEG95 spectra. They ob- 
tained encouraging results, which showed that forbidden line lumi- 
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Wavelength [Angstrom] 

Figure 6. OI 6300 A equivalent width histogram for the LVC sources in 
HEG95 compared to the model predictions for primordial discs (red line) 
irradiated by X-ray luminosities between 2xl0 28 and 2xl0 30 erg/sec, and 
compared to inner hole sources (blue line) with hole sizes between ~8 and 
~30 AU irradiated by an X-ray luminosity of 2xl0 30 erg/sec. The error in 
the equivalent widths is around the 10-20% level. 



nosities roughly of the observed magnitude can be produced in the 
ionised wind from an EUV-photoevaporated disc, and highlighted 
instead the failure of magnetic wind models in producing these 
lines. They identified, however, a discrepancy between their model 
and the observations in the luminosity of the OI 6300 A line which 
remained underpredicted by approximately two orders of magni- 
tude in their models. 

The failure of the Font et al (2004) calculations to repro- 
duce the observed OI 6300A line luminosities results from the 
fact that an EUV-driven wind is, by construction, fully ionised. 
The OI 6300A line can only then be produced in a thin layer 
at the bottom of the photoevaporating (ionised) layer. An X-ray 
driven photoevaporation process, however, yields a disc wind that 
is predominantly neutral, thus allowing abundant production of the 
OI 6300A line. Indeed we predict luminosities of a few 10 ^ L© 
for h x = 2xl0 30 erg/s, roughly two orders of magnitude higher than 
the value predicted by Font et al. (2004) and in the range of values 
reported by HEG94. Figure [6] shows a comparison of equivalent 
widths observed in the LVC of OI 6300A by HEG95 (black his- 
togram) to the those predicted by our models (red and blue line). 
The red line represents the primordial disc models irradiated by X- 
ray luminosities between 2xl0 28 and 2xl0 30 erg/sec and the blue 
line corresponds to inner hole sources with hole sizes between ~8 
and ~30 AU irradiated by an X-ray luminosity of 2xl0 30 erg/sec. 
The higher equivalent widths shown on the histogram are lines 
where even the LVC velocity range (-60 to 60 km s~') is saturated 
by emission from a jet (the x-axis is truncated in this plot and could 
go up to 7A). We note however that the range of hole sizes that 
we have investigated is restricted, furthermore we do not have any 
inner-hole models for lower luminosity. 

The OI 6300A line is produced by collisional excitation, with 
neutral hydrogen being the dominant collider (rates from Launay 
& Roueff, 1977). We stress that in T-Tauri stars with winds colli- 
sional excitation of the 'D 2 level of O will be the dominant pro- 
duction mechanism of the OI 6300A line and not OH dissociation. 



UV photodissociation of OH, which may lead to branching to the 
excited 'D2 level of neutral O, producing fluorescent emission of 
the OI 6300A and OI 6363A lines (Storzer & Hollenbach 2000) is 
generally ruled out for the HEG95 sample by the fact that too large 
FUV fields are required and because of the difficulty in reproducing 
the observed OI 6300 A / OI 5577 A ratio (see also discussion in 
Storzer & Hollenbach 2000 and Font et al. 2004). 

The luminosity of the OI 5577 A is underpredicted by our 
models compared to the HEG95 observations. This is, however, not 
a cause of concern, since this is most likely due to the fact that 
collisions with neutral hydrogen (which are likely to dominate as in 
the case of the OI 6300 A line) are not included in our calculation, 
due to the unavailability of collision rates for the 1 So level (Launay 
& Roueff, 1977, only include the first four levels, 3 P2, 3 Pi, 3 Po and 
'D 2 ). The OI 2973 A line is also a lower limit in our model for the 
same reason. 

In brief, an X-ray driven wind model offers a natural route to 
the production of the observed LVC of the OI 6300 A emission line. 
In the next section we will discuss the profiles we obtained for this 
line and compare them with observations. 

Finally, Oil and OIII forbidden lines are also listed in Tables[T] 
and [3] OIII forbidden lines are not very strong for primordial disc 
models and therefore are only listed for inner-hole models. Inner 
hole model winds are less dense and hotter than those of primordial 
discs, hence favouring the production of nebular-like lines like the 
Oil and OIII listed in the table. We draw particular attention to 
the OIII 88.3yum line, which has luminosities of around 10~ 7 L o for 
inner hole sources, and may be just detectable with PACS for the 
nearer sources. 

4.3 Neon 

Since the prediction by Glassgold et al (2007) and the first detec- 
tions of the Nell 12.8 fim line in T-Tauri discs (Pascucci et al 2007; 
Espaillat et al 2007; Lahuis et al 2007; Herczeg et al 2007), this line 
has received a lot of attention in the literature, both observationally 
(Giidel et al 2010, submitted; Pascucci & Sterzik 2009; Najita et 
al 2009; Pascucci et al 2009; van Boekel et al 2009; Flaccomio et 
al 2009) and theoretically (Meijerink et al, 2008; Schisano et al. 
2009; Alexander 2008; Hollenbach & Gorti 2009). This line was 
originally proposed as a new probe of X-ray irradiation of T-Tauri 
discs, but a clean correlation with X-ray luminosity does not seem 
to exist even in the latest data. The large scatter in the L(NeII) for 
stars of a given Lx may be due to differences in the irradiating spec- 
trum or the disc structure (Schisano et al 2009), variability, or, for 
the higher L(NeII) cases, to stellar jets (Hollenbach & Gorti 2009; 
Shang et al 2010, submitted). 

The L(NeII) we calculate from our models are in the range of 
those observed and show a roughly linear dependance with X-ray 
luminosity (see also Section [3TT ). What our models show, however, 
is that this line is produced almost exclusively in the photoevapo- 
rating wind, as clear from the blue-shifted profile for non-edge on 
inclinations (Tableland FigureQJ. This opens the exciting prospect 
of using this line as a probe of photoevaporative winds as we will 
discuss further in the next Section. It is also clear that a clean linear 
relationship with L x is not expected (or found) given that the wind 
structure and the Nell emitting region change in response to differ- 
ent X-ray luminosities. In this perspective it is easy to understand 
the scatter in the L(NeII)-Lx relation noticed by Giidel et al (2010, 
submitted) even for the non-jet sources. 

L(NeII) is about a factor of two larger for inner-hole models 
than a primordial model of the same X-ray luminosity and this is 
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simply due to the larger emission region of the inner-hole sources 
(see Figure[3j. We note that while Schisano et al. (2009) used disc- 
only models to probe the dispersion in the L(Ne)-Lx relationship, 
many of their conclusions with regards to the nature of the scatter 
(e.g. hardness of the ionising radiation field) are equally applicable 
to a line that is produced in the wind. 

The Nelll to Nell ratio for our models is never lower than 
~0. 1 , while observations generally indicate lower values, e.g. -0.06 
(Lahuis et al 2007), -0.045 (Najita et al 2010), although in the case 
of WL5 a value of -0.25 is reported by Flaccomio et al (2009). 
Unfortunately only very few detections of this line are currently 
available and it is therefore hard to draw statistically significant 
conclusions. To make matters worse this ratio is quite uncertain 
in theoretical models due to the lack of neutral collision rates, as 
mentioned in Sectionf37T] and uncertainties in the charge exchange 
coefficients. In view of the current observational and theoretical 
difficulties with this line ratio, its diagnostic power cannot yet be 
harnessed. 

4.4 Magnesium 

We predict a rather strong Mgl line at 4573 A (~10~ 5 L o for 
L x = 2xl0 30 erg/sec). This line has been detected to date, in three 
young low mass stars, namely VY Tauri (M0; Herbig 1990), XZ 
Tau B (Ml; White & Ghez 2001) and very recently in the low-mass 
variable TWA 30 (Looper et al. 2010, submitted). Detection of this 
line in M-stars is rendered possible by the much lower stellar con- 
tinuum in the B-band compared to higher mass stars. The strength 
of this line implies that it could be an important coolant and should 
be considered when thermal calculation are carried out. 

4.5 Silicon 

We include Si lines in Tables [TJ and [3] and draw the reader's at- 
tention to the Sill 34.8pm transition which has a luminosity of 
approximately 10~ 6 L o for h x = 10 30 erg/sec. We caution however 
that Si lines and Mg lines are strongly affected by depletion onto 
dust grains which is very poorly constrained. While a hindrance 
for the prediction of accurate line luminosities, the dependance of 
this lines on depletion factors may actually be used to obtain clues 
about the disc evolution through comparison of models with obser- 
vations. 

4.6 Sulfur 

We predict a number of SII, SIII and SIV lines all with fairly large 
luminosities. In particular our models predict SII 67 17A and 673 1 A 
lines in the range observed by HEG9fu The ratio of the 6731 A to 
the 6717A, which is a useful density diagnostic, is approximately 
two in our models, and it compares well with the value reported for 
most of the sources in HEG95. A ratio of two at -7000 K corre- 
sponds to an electron density of just below 10 4 cm~ 3 , again consis- 
tent with a wind origin of this line. Furthermore, all of the detec- 
tions for the LVCs of these lines in HEG95 show a small blue-shift, 
reinforcing the conclusion that all discs show evidence of a wind 
which can be produced by photoevaporation. The OI to SII ratios 
also compare favorably, when the OI line luminosites observed are 
in the range predicted by the model (i.e. < 10 _4 L o ), however we 
note that this ratio can be as large a one order of magnitude for the 

5 Note that here we only compare with HEG95 for the LVC detections 



higher OI luminosities detected by HEG95, these detections corre- 
spond to objects where the jet emission has saturated the LVC and 
are associated with large equivalent widths in the LVC OI line (e.g. 
CW Tau, DQ Tau and HN Tau). 

4.7 Hydrogen 

Tables [2] and [4] show the luminosities of the hydrogen recombi- 
nation lines. The HI(7-6) and HI(6-5) transition were detected by 
Pascucci et al. (2007), Ratzka et al. (2007), and a larger set of IR 
lines have been detected from the transition object TW Hya by Na- 
jita et al. (2010). The reported luminosities of those lines are of 
the same magnitude as for the Nell 12.8/jm line also detected in 
the same objects. Our models, on the other hand, can only pro- 
duce hydrogen recombination lines that are almost two orders of 
magnitude less luminous than the Nell 12.8yum. It seems therefore 
unlikely that a disc or wind origin can explain the reported hydro- 
gen recombination line luminosities. This had already been noticed 
by Hollenbach & Gorti (2008), who analytically derive L(HI 7- 
6)/L(NeII)~0.008, similar to what we find (0.015). We agree with 
Hollenbach & Gorti's conclusion that these IR hydrogen lines must 
originate in a region of very dense plasma (where Nell 12.8yum 
would be collisionally quenched) close in to the star, possibly in an 
accretion shock. This is in line with earlier works (e.g. Hartmann 
et al 2004; Muzerolle et al 1998, 2000; Kurosawa et al 2005) that 
H recombination lines mostly originate in the magnetospheric ac- 
cretion flow. 



5 DISCUSSION: BLUE-SHIFTED LINE PROFILES, A 
SMOKING GUN FOR PHOTOEVAPORATION 

The recent observation of a blue-shifted Nell 12.8yum line profile by 
Pascucci et al. (2009, P09) in a few well known YSOs with discs, 
including TW Hya, has been interpreted as evidence for photoevap- 
oration. As shown by Alexander (2008) the EUV-wind solution of 
Font et al. (2004) produces a good match to the data (red line in Fig- 
ure|7J, as does the X-ray driven wind of O10 discussed in this paper 
(blue line in Figure [7). The X-ray wind model used for this figure 
is that of the primordial disc model irradiated by an X-ray luminos- 
ity, Log(Lx) = 30.2 erg/sec. We stress that neither the EUV-wind 
or X-ray wind models were statistically fit to the data to produce 
the match. In both cases a model was chosen from the set of exist- 
ing ones that best reproduced the observations. Both the EUV wind 
and the X-ray wind models initial parameters could be tweaked to 
produce a perfect fit to the data, but this was not attempted in this 
paper or by Alexander (2008) as the point of Figure [7j is only to 
show that profiles in the range of those observed can be reproduced 
by photoevaporation models. 

As both the EUV- and X-ray model produce a good match 
to the Nell data, it seems that this particular diagnostic cannot be 
used to discriminate between them. The OI 6300A line intensities 
and LVC profiles, on the other hand, can only be reproduced by 
an X-ray driven photoevaporation model, which predicts a largely 
neutral wind. 

While a blue-shifted emission line profile is a clear tell-tale 
sign of outflows, the detection of non-blue-shifted lines (e.g. Na- 
jita et al. 2009) cannot be unambiguously interpreted as absence 
of a wind. It is clear from Figure Q] and [2] and the discussion in 
the previous section, that there are many reasons why a line like 
Nell 12.8yum or OI 6300A, that has a significant blue-shift in some 
cases, may instead appear centred in others -e.g. inclination, inner 
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Figure 7. Profile of the Nell 12.8,um fine-structure line as predicted by the 
X-ray driven photoevaporation model of O10 -blue line-, the EUV-only 
model of Font et al. 2004 (see Alexander 2008) -red line-, compared to the 
observation of TW-Hya by Pascucci et al. (2009) -black squares-. 



hole size, X-ray luminosity-, also different profiles respond differ- 
ently to the change in these parameters, such that the relative cen- 
troid of, say, the Nell 12.8yum with respect to the OI 6300A line 
may also shift with the shift in X-ray luminosity (see Figure 1). 

It should also be made clear that photoevaporation models pre- 
dict that a large number of primordial discs, not only so-called tran- 
sition discs, should have a disc wind and therefore produce blue- 
shifted emission lines. This is in agreement with HEG95's obser- 
vations that showed blue-shifted LVC of OI 6300A for all the YSOs 
with discs. The detection of blue-shifted lines does not indicate that 
a given disc is evolved, all that it indicates is that energetic radi- 
ation from the central star is able to reach the disc establishing 
a flow. Our models predict intensities and blueshifts for this line 
that are in agreement with observations, suggesting that the LVC of 
OI 6300A is indeed a good candidate for tracing photoevaporating 
disc gas. Further spectrally resolved observations of this line would 
be extremely helpful in testing the presented model. 

Pascucci et al. (2009) found that of the seven sources they ob- 
served with VISIR, a Nell detection was only obtained in the 'tran- 
sition' (i.e. dust inner hole) sources. They interpreted their result 
as evidence of EUV-driven photoevaporatior^ and concluded that 
this can occur only at a late stage in the disc evolution. However 
we note that the detection limits at the distance of the primordial 
disc sources which showed no evidence of Nell emission are of a 
few 10~ 6 L G , these are quite high and would only allow the strongest 
Nell sources to be detected. Indeed Giidel et al (2010, submitted) 
show that, excluding sources with known jets, many cTTS have 

6 At the time of publication of Pascucci et al. (2009) line profiles had only 
been calculated for the EUV-driven wind model (Font et al. 2004, Alexander 
2008). 



L(NeII) lower than a few 10~ s L o . The Nell luminosities predicted 
by our models for a primordial disc irradiated by Log(Lx)=30.2 
would be barely detectable at the sensitivity of the Pascucci et al. 
(2009) observations. The fact that Nell 12.8//m was only detected 
in three of the sources of the P09 sample and that these happened 
to be transition discs may just be a coincidence, or, more likely, it 
may be due to the fact that, as predicted by our models, the Nell 
line is brighter for transition discs. The latter explanation however 
would only stand if the inner hole of these sources were devoid of 
gas (or at least the gas was optically thin to both X-rays and the 
interested transitions). However all of these 'transition' discs are 
still accreting and depending on the gas densities in the inner hole, 
perhaps a primordial disc model may be more appropriate. 

One should also point out that the magnitude of the blue-shift 
depends entirely on the size and location of the emitting region. The 
blueshift will be largest for inclinations where the emitting region is 
such to maximise the line-of-sight component of the velocity vec- 
tors. The blueshift of a line is therefore not necessarily maximised 
for a completely face on disc. A good example of this is provided 
by the OI 6300A profile of the Log(L^) = 29.3 erg/sec primordial 
disc model. 

While the current observational and theoretical emission line 
profile sample is still too small to draw any significant conclusions, 
these observations highlight the potential of using high resolution 
spectroscopy of YSOs and detailed profile modelling to learn about 
the dispersal mechanism of protoplanetary discs. 



6 CONCLUSIONS 

X-rays have recently been shown to drive powerful outflows (Er- 
colano et al. 2008b, 2009a, Owen et al 2010a). The recent radiation- 
hydrodynamic calculation of Owen et al. (2010a) have demon- 
strated that that the two time-scale dispersal of protoplanetary discs 
around solar-type stars is successfully reproduced by this model. 
Further evidence is provided by the observation of shorter disc life- 
times in low-metallicity environments (Yasui et al. 2009), a be- 
haviour that is expected if X-ray driven photoevaporation dom- 
inates the dispersal of discs and argues against planet formation 
as the main dispersal mechanism (Ercolano & Clarke, 2009). Fur- 
thermore the X-ray spectra of T-Tauri stars are fairly well charac- 
terised (e.g. Giidel et al 2007) and the observed luminosities imply 
that significant X-ray photoevaporation must occur, unless the in- 
ner source is completely obscured due to accretion funnels thicker 
than 10 22 cnr 2 (Ercolano et al. 2009a) and covering all lines of sight 
from the source to the disc. 

Encouraged by these results, in this paper we have presented 
an atlas of emission lines from atomic and low-ionisation species 
emitted by photoevaporating protoplanetary discs around solar type 
stars. We have used disc and wind density distributions obtained by 
the radiation-hydrodynamic calculations of Owen et al (2010a) in 
the X-ray photoevaporation paradigm and have considered discs 
irradiated by a range of X-ray luminosities and at various stages of 
clearing, from primordial (i.e. extending all the way into the dust 
destruction radius) to inner hole sources with various hole sizes. 
Line profiles for the Nell 12.8yum, OI 6300A and HI recombination 
lines have also been calculated and should provide guidance for the 
interpretation of high-resolution spectroscopy of YSOs. 

We compare the predictions of our models to available data 
in the literature and find agreement for the line intensities and pro- 
files. We draw particular attention to the good agreement with the 
observations of Hartigan et al (1995) of the OI 6300A line intensity 
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and profile in the spectra of T-Tauri stars. This had been a problem 
for previous, EUV-driven photoevaporation models (e.g. Font et al. 
2004) which underestimated its luminosity by a couple of orders of 
magnitude. The problem with these previous models lay in the fact 
that EUV winds are, by construction, ionised, while X-ray photo- 
evaporation provides a mechanism to drive warm neutral winds, a 
condition necessary for the production of blue-shifted collisionally 
excited 01 6300A and OI 5577A. 

As noted by Font et al (2004), magnetic wind models system- 
atically underpredict the luminosities of forbidden lines and cannot 
fit the observed profiles (e.g. Garcia et al 2001). The main problem 
with magnetic wind models is that they are too cold - an X-ray 
driven photoevaporative wind naturally produces temperature, ion- 
isation degrees and consistent kinematic patterns able to produce 
forbidden line intensities and profiles in agreement with the obser- 
vations. These results suggest that, while magnetic disc winds may 
work in conjunction with photoevaporation, they are unlikely to be 
operating alone. 

We analyse our results for the Nell 12.8/mi line profiles in 
the context of recent high resolution observations of Pascucci et al 
(2009) and Najita et al. (2009). We discuss how the detection of a 
blue-shifted line is a clear sign of an outflow, but caution against 
interpreting non-blue-shifted profiles or non-detections as absence 
of a disc wind. We argue that the non-detections for the primordial 
discs in the Pascucci et al (2009) sample could be explained by in- 
sufficient sensitivity of the observations that are only able to probe 
the upper end of the Nell luminosity distribution. The non-blue- 
shifted line centroids in the Najita et al. (2009) sources are probably 
due to unfavourable (high inclination) viewing angles or to the ef- 
fect of seeing blue and redshifted emitting material through a large 
inner hole, as may be the case of GM Aur, although we caution that 
this explanation is only valid if the gas, known to be present in the 
inner (dust) hole of this object, is optically thin to X-ray radiation 
and to the line in question. 

Finally, our models cannot reproduce the observed emission 
line luminosities of IR hydrogen recombination lines, which for 
some objects are reported to be as bright as the Nell 12.8/mi (Pas- 
cucci et al 2007; Ratzka et al 2007). As suggested by Hollenbach & 
Gorti (2009), these lines are probably produced in a dense, ionised 
region that is very close to the central star. High resolution spec- 
troscopy of these lines is needed to further constrain their origin. 
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